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Abstract 



We consider a transport equation of the Boltzmann-Langevin type for non-Abelian 

0^ ' plasmas close to equilibrium to derive the spectral functions of the underlying mi- 

,^ ' croscopic fluctuations from the entropy. The correlator of the stochastic source is 

'j^. obtained from the dissipative processes in the plasma. This approach, based on 

Qh. classical transport theory, exploits the well-known link between a linearized collision 

(— I \ integral, the entropy and the spectral functions. Applied to the ultra-soft modes of 

a hot non-Abelian (classical or quantum) plasma, the resulting spectral functions 



agree with earlier findings obtained from the microscopic theory. As a by-product, it 
follows that Bodeker's effective theory is consistent with the fluctuation-dissipation 



5t , theorem. 
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It has been recognized that dynamical properties of (non-perturbative) quasi-particle 
excitations in non-Abehan plasmas can be described very efficiently by means of effective 
transport equations. A prominent recent example is given by Bodeker's effective classical 
theory for the ultra-soft modes in a hot non-Abelian plasma close to equilibrium IQ, which 
corresponds to a transport equation of the Boltzmann-Langevin type. In |0, a general 
procedure has been presented, based on classical coloured point particles, to obtain effective 
transport equations from the microscopic theory after integrating-out the fluctuations about 
the mean values, and taking the Gibbs ensemble average in phase space. On the one-loop 
level, the same mean field equations of [^ have been obtained recently within a many-particle 
world line formalism Q] (see also 0]). The collision integral and the source of stochastic 
noise of [Q have been obtained from to leading order in a weak coupling expansion, and 
at logarithmic accuracy. It was also realized that the dynamical equations are the same for 
the classical and the quantum plasma, changing only in the value of the Debye mass [0]. 
Other approaches to obtain the collision term of |^ have been reported as well [§-0|. 

In the present article we consider, based on classical transport theory, a generic 
Boltzmann-Langevin equation for the one-particle distribution function f{x,p,Q), given 
as 

P' [-i^-9r'"'KQ'^a-9QaF;;,^) f{x,p,Q) = C[f]{x,p,Q) + ax,p,Q) . (1) 

Here, the variables Q describe the non-Abelian colour charges. The transport equation 
contains an effective collision term C[f] and an associated source for stochastic noise (. The 
SU(A^) gauge fields appearing in the above equation are self-consistent, that is, generated 
by the same particles of the plasma. The Yang-Mills equation are 

{D,Fna = J:i^)^9 E JdPdQQ,p''f{x,p,Q), (2) 

hclicitics 
species 

where the momemtum measure reads dP = d'^p2Q{po)5{p'^ — m?), and the colour measure 
dQ was defined in [0. The current in (H) is covariantly conserved. We work in natural 
units c = h = ks = ^, unless otherwise specified. From now on we will omit the sum over 
different species of particles and helicities. 

In the collisionless limit C = ( = 0, the above set of transport equation reduces to those 
introduced by Heinz [Q. In the general case however, the r.h.s. of (|I]) does not vanish due 
to effective interactions (collisions) in the plasma, resulting in the term C[f]. In writing 
(HI), we have already made the assumption that the one-particle distribution function / is a 
fluctuating quantity. This is quite natural having in mind that / describes a 'coarse-grained' 
microscopic distribution function for coloured point particles, and justifies the presence of 
the stochastic source ( in the transport equation. For non-charged particles, a similar 
kinetic equation has already been considered in |Q] (see also |TD|, where stochastic noise is 
introduced to a Schwinger-Dyson approach). 

Given the stochastic dynamical equation (|I]), the question raises as to what can be said 
on general grounds about the spectral functions of / and (. Here, we shall assume that 
the dissipative processes are known close to equilibrium, but no further information is given 



regarding the underlying fluctuations. This way of proceeding is complementary to [Q, where 
the r.h.s. of (0) has been obtained from correlators of the microscopic statistical fluctuations. 
We then show that the spectral function of the fluctuations and the noise correlator close to 
equilibrium can be obtained from the knowledge of the entropy of the plasma, and from the 
dissipative term in the effective transport equation. This gives a well-deflned prescription 
as to how the correct source for noise can be identifled without the detailed knowledge 
of the underlying microscopic dynamics responsible for the dissipation. The basic idea 
behind this approach relies on the essence of the fluctuation-dissipation theorem (FDT). 
While this theorem is more general, here we will only discuss the close to equilibrium 
situations. According to the FDT if a fluctuating system remains close to equilibrium, then 
the dissipative process occurring in it are known. Vice versa, if one knows the dissipative 
process in the system, one can describe the fluctuations without an explicit knowledge of 
the microscopic structure or processes in the system. The cornerstone of our approach is 
the entropy of the fluctuating system, which serves to identify the thermodynamical forces, 
and leads to the spectral function for the deviations from the non-interacting equilibrium. 
Before entering into the discussion of plasmas, we will illustrate this way of proceeding by 



reviewing the simplest setting of classical linear dissipative systems |jTT|] . A generalization to 
the more complex case of non-Abelian plasmas will then become a natural step to perform. 
We consider a classical homogeneous system described by a set of variables Xj, where i is 
a discrete index running from 1 to n. These variables are normalized in such a way that 
their mean values at equilibrium vanish. The entropy of the system is a function of the 
quantities Xi, S{xi). If the system is at equilibrium, the entropy reaches its maximum, and 
thus {dS/dxi)eq = 0, Vz. If the system is taken slightly away from equilibrium, then one can 
expand the difference AS = S — Seq, where S^q is the entropy at equilibrium, in powers of 
Xi. If we expand up to quadratic order, then 

The matrix (3ij is symmetric and positive-definite, since the entropy reaches a maximum at 
equilibrium. The thermodynamic forces Fi are defined as the gradients of AS 

For a system close to equilibrium the thermodynamic forces are linear functions of Xi, 
Fi = PijX^ . If the system is at equilibrium, the thermodynamic forces vanish. In more 
general situations the variables Xi will evolve in time. The time evolution of these variables 
is given as functions of the thermodynamical forces. In a close to equilibrium case one can 
expect that the evolution is linear in the forces 



dx^ 
It 
which, in turn, can be expressed as 

(ix* 
Itt 



-l''F, + C , (5) 



-X'^x, + C , (6) 



The first term in the r.h.s. of the above equation describes the mean regression of the system 
towards equilibrium, while the second term is the source for stochastic noise. The quantities 
Y^ are known as the kinetic coefficients, and it is not difficult to check that •jij = XikPL^- 
From the value of the coefficients jSij one can deduce the equal time correlator 

{x,it)x,it)) = (3-^ , (7) 

which is used to obtain Einstein's law 

(x\t)F,it))=S]. (8) 

After taking the time derivative of (H), assuming that the noise is white and Gaussian 

(Cit)C^{t')) = u^^it-t'), (9) 

we find that the strength of the noise self-correlator u is determined by the dissipative 
process 

z/^j = yi + y^, (10) 

which is the FDT relation we have been aiming at. 

We now come back to the case of a non-Abelian plasma and generalize the above dis- 
cussion to the case of our concern. We will consider the non-Abelian plasma as a linear 
dissipative system, assuming that we know the collision term in the transport equation. 
In order to adopt the previous reasoning, we have to identify the dissipative term in the 
transport equation, and to express it as a function of the thermodynamical force obtained 
from the entropy. The deviation from the equilibrium distribution is given here by 

A/(x,p,g) = /(x,p,g)-/eqM, (11) 

and replaces the variables Xi discussed above. The entropy flux density for classical plasmas 
is given as 

S,{x) = - J dPdQp,f{x,p,Q) {\n{f{x,p,Q)h') - l) , (12) 

where h is an arbitrary constant such that f{x,p, Q)h^ is dimensionless. The /z = compo- 
nent of (p!2| ) gives the entropy density of the system. The entropy itself is then obtained as 
S = Jd^xSo{x). 

We shall now assume that the deviation of the mean particle number from the equilibrium 
one is small within a coarse-graining volume. This can always be arranged for at small 
gauge coupling (7 ^ l.[| We then obtain AS' just by expanding the expression of the entropy 



••■More precisely, in the close-to-equilibrium plasma we consider / as being coarse-grained over a 
Debye volume. This entails that (i) the particle number fluctuations within a Debye volume are 
parametrically suppressed by powers of the plasma parameter (that is, ultimately, by g), (ii) the 
two- and higher-particle distribution functions are suppressed as opposed to /, and {Hi) typical 
relaxation scales or times are not affected by the coarse-graining. 



density in powers of A/ up to quadratic order. It is important to take into account that we 
will consider situations where the small deviations from equilibrium are such that both the 
particle number and the energy flux remain constant, thus 

JdPdQ<^{p)Af{x,p,Q) = 0, for <l>(p) = po, PoPm • (13) 

Under those assumptions, one reaches to 

ASo(x) = - / dPdQpo ^j. , . 

J 2/cq(u;p) 

where in the last equality we have taken into account the mass-shell condition, with po = 
LOy = a/p^ + m'^. Without loss of generality, we will consider from now on the case of 
massless particles, so Up = p = |p|. 

The thermodynamic force associated to A/ is defined from the entropy as 

SAf{x, P, Q) /eq(p) 

We now linearize the transport equation (|ip and express the collision integral close to equilib- 
rium in terms of the thermodynamical force. Dividing (|lD by po and imposing the mass-shell 
constraint, we find 



'^o dp 



v^D.Af - gv^Q^F:,^ = C[Af]ix, p, Q) + C(x, p, Q) , (16) 



where v^ = p'^/po = (l,v), with v^ = 1. We also introduced the shorthand D^Af = 
{dfj, — gf"'^'^A^^bQcd^)Af 0. It is understood that the collision integral has been linearized, 
and we write it as 

C[Af] (t, X, p,Q)=J d'x'd'p' dQ' K(x, p, Q; x', p', Q')Af{t, x', p', Q') , (17) 

with t = Xq. For simplicity, we take the collision integral local in time, but unrestricted 
elsewise.| According to the FDT, the source of stochastic noise has to obey 

(c(.. p. «c(.'. P'. Q')) ^ - i'-§PP^ + Wr^ ) (i«) 

V 6F{x',p',Q') 6F{x,p,Q) J 



in full analogy to (10). With the knowledge of the thermodynamical force (15) and the 



linearized collision term (\[7\) we arrive at 



^Of course, gauge invariance imposes further conditions on both the collision term and the noise. 
However, these constraints are of no relevance for the present discussion. 



iCix, p, Q)C(x', p', Q')) = - (/eq(p)i^(x, p, Q; x', p', Q') + sym.) 6(1 - f) . (19) 

Here, symmetrisation in (x, p, Q) <-> (x', p', Q') is understood. 

Notice that we can derive the equal time correlator for the deviations from the equilib- 
rium distribution simply from the knowledge of the entropy and the thermodynamical force, 
exploiting Einstein's law in full analogy to the corresponding relation (|^). Using (^) we 
find 

(A/(x, p, Q) A/(x', p', Q'))*=i' = h^{v)5^'^ (x - x0(5(3) (p - p')<5(Q - Q') . (20) 



If the fluctuations A/ have vanishing mean value, then (^0]) reproduces the well-known result 
that the correlator of fluctuations at equilibrium is given by the equilibrium distribution 
function. In order to make contact with the results of 0], we go a step further and consider 
the case where A/ has a non-vanishing mean value to leading order in the gauge coupling. 
Splitting A/ = gf^^^ + 5f into a deviation of the mean part (A/) = gf^^^ and a fluctuating 
part {5f) = and using (EUD, we obtain the equal time correlator for the fluctuations Sf as 



{6 fix, p, Q)6fix', p', Q'))*=*' = fe^ip)S^'\^ - x')5(3)(p - p')6{Q - Q') 

-//«(x,p,g)/«(x',p',Q'). . ■ (21) 



t=t' 



This result agrees with the correlator obtained in ||^ from the Gibbs ensemble average as 
defined in phase space in the limit where two-particle correlations are small and given by 
products of one-particle correlators. 

Up to now we have dealt with purely classical plasmas. On the same footing, we can 
consider the soft and ultra-soft modes in a hot quantum plasma. These can be treated 
classically as their occupation numbers are large. The sole effect from their quantum nature 
reduces to the different statistics, Bose-Einstein or Fermi-Dirac as opposed to Maxwell- 
Boltzmann. The corresponding quantum FDT reduces to an effective classical one |11,12|. 

Some few changes are necessary to study hot quantum plasmas. As in [^, we change 
the normalisation of / by a factor of (27r/i)'^ to obtain the standard normalisation for the 
(dimensionless) quantum distribution function. Thus, the momentum measure is also mod- 
ified by the same factor, dP = d'^p2Q{po)6{p'^) / {27ihY for massless particles, and ^ = 1. To 
check the FDT relation in this case one needs to start with the correct expression for the 
entropy for a quantum plasma. The entropy flux density, as a function of f{x,p, Q), is given 
by 

S,{x) = - j dPdQp^ (/In/ T (1 ± /) In (1 ± /)) , (22) 

where the upper/lower sign applies for bosons/fermions. From the above expression of the 
entropy one can compute AS*, and proceed exactly as in the classical case, expanding the 
entropy up to quadratic order in the deviations from equilibrium. Thus, we obtain the noise 
correlator 

(C(x, p, Q)C(x', p', Q')) = ~{2nf (/eq(p)(l ± /cq(p))i^(x, p, Q; x', p', Q') + sym.) 5{t - t') . 

(23) 



Again, the spectral functions of the deviations from equihbrium are directly deduced from 
the entropy. As a result, we find 

{Af{x,p,Q)Af{x',p\Q%^, = {2nff,^{p){l±Up))6^'\^-^')6^'\p-p')6{Q-Q') . 

(24) 

Expanding A/ = gf^^^ + Sf as above, we obtain the equal time correlator for Sf, which 
agrees with the findings of in the case where two-particle distribution functions can be 
expressed as products of one-particle distributions. 

With the knowledge of the above spectral functions for the fluctuations in a classical or 
quantum plasma one can derive further spectral distributions for different physical quanti- 
ties. In particular, we can find the correlations of the self-consistent gauge field fluctuations 
once the basic correlators as given above are known. This is how those spectral functions 
were deduced in |Q]. 

As a particular example of the above we consider Bodeker's effective kinetic equations 
which couple to the ultra-soft gauge field modes with spatial extensions ^ l/rrijj, where 
irto is the Debye mass in a non-Abelian plasma close to equilibrium. The linearized collision 
integral has been obtained to leading logarithmic accuracy by several different approaches 
|]1|J^J^-|^ . To leading order, they all employ an IR cut-off of the order of gm£, for the elsewise 
unscreened magnetic sector. 

We will first concentrate on the classical plasma, for particles carrying two helicities. 
It is most efficient to write the transport equation not in terms of the full one-particle 
distribution function, but in terms of the current density 

J,'ix,v)=Angv' JdpdQp'QaAfix,p,Q). (25) 

(Notice that /eq gives no contribution to the current.) The current of (Q) follows after 
integrating over the angles of v, J^{x) = J ^J^a{x, v) 0. Expressed in terms of (pSj), the 
linearized Boltzmann-Langevin equation ([l6|) becomes 

KD^, J']{x, v) = -mlv'v^F^%x) + v''C[j'']{x, v) + C{x, v) . (26) 

where tud is the Debye mass [Q 

ml = -8ng'C,ldpp'^ (27) 

and the constant C2 depends on the representation of the coloured particles 

j dQQaQb = C25ab. (28) 

The linearized collision integral is related to ( |T7|) by 

C[J^] {x, w)=ATig j d^x' dfiv' dp dp' dQ dQ'p'p'^ Qa K{x, p, Q; x', p', Q') A/(t, x', p', Q') 

(29) 



and has been obtained explicitly [l|,0,|]-0] as 

dQv' 



C[X]i^,^r) = -7 /^X(v,v')J-i'(x,v') (30) 



where the kernel reads 

J(v, v') = 5(2) (V - v') - 1 (I^l2^ (31) 

^ .^1 - (V ■ V')2 

and 7 = g"^ NT In {l/g)/A7[. Comparing (p9|) with (^) we learn that only the part of the 
kernel K which is symmetric under (x, p, Q) ^->- (x', p', Q') contributes in the present case. 
This part can be expressed as 

K(x, p, Q; x', p', g') = -7 ^^^Sip - p')5{Q - Q')5^'\^ - x') . (32) 

According to our findings above, the self-correlator of the stochastic source for the classical 
plasma obeys 

(C(a;, v) a{y, v')) = (47r)2 g^ J dpdp' dQ dQ' p'p" Q^Q', v^v"" (C(a;, p, Q) ((y, p', Q')) 

= 2 7 T m^ v'^v"' J( V, v') d^b 5^^ {x ~ y) . (33) 

The sum over the two helicities of the particles has been taken into account. In order to 
obtain (0), we have made use of (|19D, (^) to (^Of), and of the relation /eq = —Tdf^q/dp 
for the Maxwell-Boltzmann distribution. 

The quantum plasma can be treated in exactly the same way. To confirm (^), we only 
need to take into account the change of normalization as commented above, and the relation 
/eq(l ±/eq) = —Tdf^q/dp for the Bose-Einstcin and Fermi-Dirac distributions, respectively. 

Using the explicit expression for the collision integral and the stochastic noise it is 
possible to confirm the covariant conservation of the current, DJ = |[l|J^. 

We thus found that the correlator (^3]) is in full agreement with the result of |l|,|^ for 
both the classical or the quantum plasma. While this correlator has been obtained in 
HIQI from the corresponding microscopic theory, here, it follows solely from the FDT. This 
way, it is established that the effective Boltzmann-Langevin equation found in |]I[ is indeed 
fully consistent with the fiuctuation-dissipation theorem. More generally, the important 
observation is that the spectral functions as derived here from the entropy and the FDT do 
agree with those obtained in p| from a microscopic phase space average. This guarantees 
that the formalism of [0] is consistent with the FDT. 

In the above discussion we have considered the stochastic noise as Gaussian and Marko- 
vian. This is due to the fact that the small-scale fiuctuations (those within a coarse-graining 
volume) are to leading order well seperated from the typical relaxation scales in the plasma. 
Within the microscpoic approach in p|, these characteristics can be understood ultimately 
as a consequence of an expansion in a small plasma parameter (or a small gauge cou- 
pling). More precisely, the noise follows to be Gaussian when higher order correlators be- 
yond quadratic ones can be neglected. The Markovian character of the noise follows because 



the ultra-soft modes are well separated from the soft ones, and suppressed in the collision 
integral to leading order. This way, the collision term and the correlator of stochastic noise 
are both local in x-space. Going beyond the logarithmic approximation, we expect from 
the explicit computation in that the coupling of the soft and the ultra-soft modes makes 
the collision term non-local in coordinate space. This non-trivial memory kernel should also 
result in a non-Markovian, but still Gaussian, source for stochastic noise. 

The present line of reasoning can in principle be extended to other approaches. Using 
the phenomenological derivation of (^) from 0, the same arguments as above justify the 
presence of a noise source with (^) in the corresponding Boltzmann equation |^|-0]. It 
might also be fruitful to follow a similar line based on the entropy within a quantum field 
theoretical language. An interesting proposal to include self-consistently the noise within a 
Schwinger-Dyson approach has been made recently in ||T0|. Along these lines, it might be 
feasible to derive the source for stochastic noise directly from the quantum field theory [^]. 

While we have concentrated the discussion on plasmas close to equilibrium, it is known 
that a fluctuation-dissipation theorem can be formulated as well for (stationary and stable) 
systems out-of-equilibrium [0. It can also be extended to take non-linear effects into 
account [|13|. Both the out of equilibrium situations and non-linear effects can be treated, 
in principle, with the general formalism presented in p|. 

CM. thanks the Institute for Nuclear Theory at the University of Washington for its 
hospitality and the Department of Energy for partial support during the completion of this 
work. 
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